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Abstract 



We describe tlie classical and quantum two dimensional nonlinear dynam- 
ics of large blue-detuned evanescent-wave guiding cold atoms in hollow fiber. 
We show that chaotic dynamics exists for classic dynamics, when the intensity 
of the beam is periodically modulated. The two dimensional distributions of 
atoms in (x, y) plane are simulated. We show that the atoms will accumulate 
on several annular regions when the system enters a regime of global chaos. 
Our simulation shows that, when the atomic flux is very small, a similar dis- 
tribution will be obtained if we detect the atomic distribution once each the 
modulation period and integrate the signals. 

For quantum dynamics, quantum collapses and revivals appear . For 
periodically modulated optical potential, the variance of atomic position will 
be suppressed compared to the no modulation case. The atomic angular 
momentum will influence the evolution of wave function in two dimensional 
quantum system of hollow fiber. 
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I. INTRODUCTION 



Evanescent light offers a very effective means of atom manipulation because of the strong 
induced dipole interaction between atoms and evanescent field when the field is far from res- 
onance. Blue-detuned evanescent light, that is light tuned above the atomic resonance, can 
be used to guide atoms in a hollow fiber [T]], 0, With blue-detuning, atoms are repelled 
from the high intensity field region near the fiber wall. The intensity in the evanescent 
field is significant over a distance of approximately ~ A in the hollow region. If the inner 
diameter is several microns, the motion of atoms will be influenced by the evanescent field 
over a large transverse area inside the fiber. This system is very suited to studying nonlinear 
quantum and classical dynamics with two degrees of freedom in the transverse plane of the 
fiber. The maximum confined transverse atomic velocities are only several centimeters per 
second, therefore the atomic dynamics will be governed by the Schrodinger equation with 
an electric dipole coupling to the evanescent field. There are few experimental studies of 
quantum nonlinear dynamics for systems with two degrees of freedom and most of our under- 
standing of quantum nonlinear dynamics and quantum chaos is based on simpler one degree 
of freedom systems, such as an atom in a modulated standing wave 0], or the microwave 
ionisation of hydrogen ||^. A notable example of experimental study of quantum chaos in 
two degree of freedom are the recent experiments on chaotic mesoscopc billiards It 
is clear that quite new phenomenon can appear in two degree of freedom systems p| that 
are absent in one degree of freedom systems. Given the success of atom optics in providing 
tests of quantum chaos in one dimension it is advisable to consider what might be achieved 
using similar techniques for two dimensional systems. In this paper we will consider the two 
dimensional quantum and classical transverse motion of cold atoms in a hollow fiber with a 
periodically modulated evanescent field. 



2 



II. THE OPTICAL POTENTIAL AND HAMILTONIAN 



A two-level atom interacting with far-off-resonant inhomogeneous laser field has an ef- 



fective potential |Tn] of the form 



U{r) = ^ln{l+p), (1) 

where A is the detuning and p = ^2^/^2/4 is a saturation parameter, with the Rabi frequency 
Q. For far-off-resonant donut beams, p <ti I, and thus 

Uir) = (2) 

In experiments |l|], 0, this potential has been used to guide cold atoms by reflecting them 
from evanescent light fields on the surface of the glass. For light striking the glass-vacuum 
interface at angle 6, the evanescent intensity profile is 

/(r) = /(O)a^ exp[2K(r - n)], (3) 

where /(O) is the input laser intensity at the fiber entrance, ri is the inner radius of fiber , 
and the factors a and k are given in terms of the index of refraction n, inner reflection angle 



9 and laser wavelength A by a = 2^n^/{rfi — l)cos6', and k = (27r/A) y (n^ sin^ 6' — 1). The 
Hamiltonian in transverse (x, y) plane for the system is 

2 2 

Ho = + K exp[2^{r - n)], (4) 

where K = ^^^Q^^ and Is is the saturation intensity. 

When the laser intensity is periodically modulated, /(O) becomes time dependent, 
/(0)[1 + ecos(co't)] and the dynamics of the atom can be chaotic for certain initial con- 
ditions. There are two factors to consider in treating the two dimensional classic chaotic 
dynamics for evanescent wave guided atoms in a hollow fiber. Firstly, the boundary of the 
hollow fiber will limit the divergence of the atomic motion in the radial direction making it 
is easier to simulate the two-dimensional the dynamics of the atomic system (and also easier 
to probe the two dimensional distribution in experiment). Secondly, due to the small inner 



radius of a hollow fiber, atoms with a small transverse velocity will be selected at the entry 
to the hollow fiber and it is not necessary to pre-cool the transverse temperature. 

The effective optical potential depends on the inner radius ri and decay coefficient n. If 
the ri is much larger than wavelength A or ^ is too large, the potential will be very steep 
on the boundary and decay rapidly away from the inner surface of the fiber. In order to 
observe two dimensional chaotic dynamics in {x, y) plane, ri should be the same magnitude 
as so that the potential will decay slowly away from the surface of the hollow fiber. 
It is possible to make the reflection angle approach the critical reflection angle sm~^(l/n) 
by using several micro hollow fibers to guide the atoms and polish the incoming end of the 
hollow fiber at a specific acute angle [Q]. 

We take ri = 2/im and 6 = 45° and consider helium as an example. The parameters 
for helium are, line width r/27r = 1.6MHz, mass M = ArUp, wave length A = 1.083/im, 
saturation intensity Ig = Q.lQmW/cim? and recoil velocity Vr = 9cm/ s. We now define 
dimensionless parameters {x,y) = {2K,x,2K,y), r = 2Kr, {px,Py) = (IfS^'l?^^)' H = 
Hj^, u) = uj/uq and i = uo^t, where ujq is a reference frequency. Omitting the tildes and 
defining ^ = ^^jl , the Hamiltonian can be rewritten as 

if(t) = ^2-±ii + eev/^-'^i(l + ecoscut), (5) 
with the canonical communication relations 

[qj,Pk\ = ikSjk, (6) 

where qj,Pk represents x, y and k = ^|^, plays the role of a dimensionless Planck constant. 
If the dimensionless Plank constant k approaches 1 , the atomic motion in small hollow fibers 
will be correctly described by quantum dynamics. 

III. CLASSIC CHAOTIC DYNAMICS AND 2D DISTRIBUTION 

Although the flux of guided atoms entering the fiber is very small, we still have the 
possibility to observe classic chaotic dynamics by integrating the signal over many periods 
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of the modulation (fig. 1). The integration length can be taken long enough to ensure that 
it includes large numbers of atoms. We will discuss the spatial distribution of the atom 
when the classical motion is chaotic. Using Hamilton's equations we find that the motion 
in the transverse plane without modulation (e = 0) is described by the equations, 

P. = -e , ev^-s (7) 

va; + y 



Py = -^y^eV^-r^ (8) 
X = Px, (9) 

y = Py (10) 

Clearly there is only one stable fixed point on the axis (r = 0). 

The choice of modulation frequency u depends on the frequency of unperturbed periodic 
motion. For simplicity we assume y = and Py = 0, so the expression for H simplifies to a 
one dimensional Hamiltonian system. The period of motion for the unperturbed Hamiltonian 



i^o is m 



I dx r^M dx , . 



where xm is determined by Hq = ^ev ^a/ '"i. Therefore 

'^O = 7^ • (12) 

JXAmo - ^e^-^^-]}-y^dx 

The graph of c^o versus Hq and xm versus Hq is shown in fig. 2. We can select the modulation 
frequency uj to control the position of the fixed points. Here we take ujq = 265KHz , 
dimensionless modulation frequency u = 2 and set ^ = 50. The new fixed points will appear 
in X ^ 0.41 and x ^ 2.87. 

We use a symplectic integration routine |T2[ , [|r^ to solve the equations of motion so as to 
preserve the Poisson bracket relation {x{t),Px{t)} = 1, and thus maintain the Hamiltonian 



character of the motion. We plot the stroboscopic portrait of the system at multiples of the 
period of modulation, t = {2n/u)s, where s is an integer referred to as the strobe number. 
From fig. 3 we can see that regions of globally chaotic motion will arise when e is large, 
together with some regular regions. A broad initial phase space distribution of atoms will 
enable some atoms to become trapped in these stable regions. 

Laser cooling and trapping techniques have the ability to cool the atom to very low 
velocities and trap them with well localized momentum, however the position distribution 
is not so well localized. Therefore the appropriate description of the initial conditions is in 
terms of a probability density on phase space {x,y,Px,Py). We define a classical state to be 
a probability measure on phase space of the form Q{x, y , p^, Py)dxdydpxdpy. The probability 
density satisfies the Liouville equation 

f = (13) 

where{,}^. p. is the Poisson bracket. This equation can be solved by the method of char- 
acteristics. To simulate an experiment, we assume atoms are initially randomly uniformly 
distributed on + y^ < Ti^, 

The momentum distributions for p^ and Py are assumed to be Gaussian distributions. 
Therefore 

Qo{x,y,Px,Py) = Qo(x)Qo(y)Qo(Px)Qo(Py), (14) 

where 

QoiPi) = exp [-{pi - p,(0)) V2apJ. (15) 

The variances of p^ and Py are related to the temperature Tj 

a,,^kBT,/[Mu;l/{2Kn (16) 

We simulated the atomic system of 10^ numbers and take = 0.1, which corresponds 

to radial rms velocity 2cm/ s for helium if ^ = 45° and u = 2 . The variance of Q function 
with time is 
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Q(r, p, t) = Qo[r(r, p, -t), p(r, p, -t)] (17) 

In the case of no modulation, atoms will accumulate around the fix point x = y = 0. 
When the modulation is added the atoms will diffuse in regions of chaotic motion but some 
will accumulate around several rings corresponding to fixed points at non-zero radius (fig. 
4). 

Because the inner size of hollow fiber is very small. The fiux of guiding atoms entering 
the fiber will be low. But we can use a large blue-detuned shield laser, with the duration 
period of square wave equals to the modulation period T = 2^/00, and the time interval 
dT << T. When atoms emit out from the hollow fiber, they will travel freely without 
the optical potential. If the shield laser turn on, it will block atoms enter the detection 
region. In the short time interval dT, atoms will enter the detection region and the detector 
records the atomic distribution of momenta and positions. This procedure can be repeated 
as long as possible till we have integrated enough atomic numbers for many snap shots. The 
simulation shows that for integration of atomic numbers at different integer strobe numbers 
, it will produce similar results (fig. 5). 

Because rings (fig. 4)represent new fixed points in phase space and the atomic radial 
momenta are almost approach zero(fig. 3). When atoms emit out and travel freely the 
shape of rings will not change very much in detection region. But the atomic velocity 
orientations outside the fixed points are randomly distributed. In order to measure the rings 
and small spatial distribution of atoms, high precision position measurements are required. 
The Raman-induced resonance imaging method El, HH, H can be used to measure the 



positions with nanometer spatial resolution limited by the quantum uncertainty principle. 

IV. TWO DIMENSIONAL QUANTUM NONLINEAR DYNAMICS 

If the transverse temperature of atoms in hollow fiber is very cold, quantum nonlinear 
dynamics will result. Here we take the same parameters and same definition for dimension- 
less parameters as the classical dynamics and the dimensionless Plank constant — 1- 



The dimensionless 2D Schrodinger equation for atoms in hollow fiber 



where 



tk^^l^ = Hm{v,t), (18) 



and 



V{x,y,t) = ^e^/^-'''{l + ecosut), (20) 



where x + y < ri . 

For simplicity we assume the boundary condition for eq. (18) is 



^(^,2/,^)l^/^5^=.,=0. (21) 



The Split Operator Method |T5|and 2D FFT (Fast Fourier Transformation) [|T^ will be 
used to obtain the numerical solution of Schrodinger eq. 18. The scheme in which the kinetic 
operator and potential operator are used to propagate the wave function separately: 

exp[-imt/k] ~ exp[-i{P)^5t/Ak] exp[-i(V)6t/2k] exp[-i(P)^5t/Ak], (22) 

and the computing errors are of 0{5t^). 

We assume at t = the wave function is a minimum uncertainty wave function. The 
initial variance of x,y are the same cr^ = ffy = u- The expression for this wave function is 

ij{x,y) = _^e-^e-^e^^«-/^e^^o«^/^. (23) 
v27r(T 

In order to observe genuine quantum nonlinear behavior we have to wait for a time that 
is longer than the classical period. The time scale for quantum nonlinear dynamics Trev is 
given by expression ^ 

Tre. = To{^\duJo/dE\y\ (24) 



In fig. 6 we have plotted the momentum mean < px > and its variance as a function of 
time for e = 0. Quantum collapses and revivals appear for quantum nonlinear dynamics as 
expected. 

In order to understand the influence of the modulation of the potential, we have plotted 
the variances of < > — < x < p^'^ > — < p^ >^ versus strobe numbers. In 
fig. 7 we found at integer strobe numbers, the modulation will increase the variance of 
momentum , but suppress the variance of position. The fluctuation of variances of positions 
and momentum are consistent with the Heisenberg uncertainty relationship in this quantum 
system. 

If we rewrite H in polar coordinator (r, 6), we get 

^ = -2(a^ + -^ + ^a^) + n^-^)^ (25) 



The general solution of eq.(19) will be K 



ij{r,e,t)=Y.a^{r,t)e'"''f', (26) 

m 

where m is angular momentum quantum number. It shows that the evolution of wave 
function will be related to atomic angular momentum. In general for t > the probability 
\ilj{x,y)f will not be symmetric except the initial wave function has no angular momentum. 
In fig. 8 for initial conditions Xq = yo = and po^ = poy = we have plotted th the 
probability at y = plane as a function of x at the time of strobe number 50, it shows that 
it will stay symmetric both for e = and 0.7. In fig. 9, for initial condition Xq = y^ = 
and pox = Poy = 1 , the probability distribution for t > will not be symmetric because 
initial atomic wave functions include angular momenta, and the interaction between atomic 
momenta and optical potential will destroy the spatial symmetry of probability. Therefore 
the atomic angular momentum plays an important rule in the evolution of wave function. 
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V. CONCLUSION AND DISCUSSION 



The atomic motion in small diameter hollow fiber is quantum mechanical if dimensionless 
Plank constant k approach 1 . Quantum collapses and revivals will appear in the dynamics 
of the mean values . If the intensity is periodically modulated, we found the modulation will 
increase the variance of momentum , but suppress the variance of position at integer strobe 
numbers. The fiuctuation of variances of positions and momentum is consistent with the 
Heisenberg uncertainty relationship in this system. In this two degree of freedom quantum 
system, the atomic angular momentum will influence the evolution of wave function. 

Although the flux of guiding atoms entering the flber is low, it's still possible to ob- 
serve the classic chaotic dynamics by integrating the signals at different strobe numbers in 
experiment. The integration length can be taken long enough to ensure that it includes a 
large numbers of atoms. We have shown that an atom moving in an intensity modulated 
evanescent-wave fleld in hollow flber can exhibit chaotic dynamics in the transverse plane. 
For atomic momenta p^, Py with Gaussian distributions, some atoms will become trapped 
in rings corresponding to radial flxed points of the modulated system in the moments of 
integer strobe numbers. 

For the atomic average radial velocity is around recoil velocity and the spatial distri- 
bution size no larger than inner diameter of flber, high resolution position and velocity 
measurements must be used. The Raman-induced resonance imaging method can reach the 
nanometer spatial resolution limited by the uncertainty principle. For atomic momentum 
distribution can be measured by Time of Flight absorption imaging method H, which 
is used to measure the temperature of super cold atoms in BEG experiment, or atomic 
velocity selection method using stimulated Raman transitions |2l|. 

It is quite interesting that both classical chaotic dynamics and quantum dynamics have 
the possibility to be realized in experiment for atoms propagating in hollow flber. Because 
the flber's inner size is very small, it will select atoms with very small radial velocity to 
enter the flber and there is no need to pre-cool the radial motion. In addition the flber's 
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boundary will limit the divergence of radial motion, making it is easier to observe both two 
dimensional classical and quantum dynamics in experiment. We believe this is a far more 
practical scheme for observing the classical and quantum nonlinear dynamics of radially 
confined atoms than other schemes such as those that use a donut beam scheme . 
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FIGURES 

FIG. 1. The diagram of proposed experiment. The large blue-detuned shield laser will block 

atoms enter the detection region except in the short time interval dT « T. The integration of 
many snap shots will have the similar effect of atomic distribution compared to one snap shot of 
large atomic numbers. 

FIG. 2. a. The relations between frequency of motion and Hamiltonian ~ Hq and xm ^ Hq. 

FIG. 3. Stroboscopic portrait of the system with e = 0.7, Pa;(0) = 0, = and y = 0. The 
maximum strobe number is 500. 

FIG. 4. The atomic distribution in {x, y) plane at the strobe number 50 for e = 0.7. The 10^ 
atoms were taken in phase space. The atoms were initially distributed on x^ -\-y^ < r\ region. The 
momenta of Px, Py are Gaussian distributions and Up^ = Gp^ = 0.1.. 

FIG. 5. The atomic distributions o/ 10^ numbers in {x,y) plane. After strobe number 50, we 
count 200 atoms at the time of every integer strobe number until s = 150. 

FIG. 6. The average momentum < Px > and momentum variance < Px > —< Px evolu- 
tions with time. The initial wave function is the least uncertainty state, ffx = cry = 0.1. For (a), 
xo = yo = 0, pox = Poy = 1-0. For (b), xq = yo = 0.5, pox = Poy = 0.0. 

FIG. 7. The position variance< x^ > — < x >^ and momentum variance < Px > —<Px>^ 
versus strobe numbers. The initial wave function is the least uncertainty state, ax = (Ty = 0.1. For 
(a) and (h) xq = yo = 0, pox = poy = 1-0. For (c) and (d), xq = yo = 0.5, pox = Poy = 0.0. 

FIG. 8. The probability distributions |^(x,0)|^ for y = 0. The initial wave function is the least 
uncertainty state, (Tx = CTy = 0.1, xq = yo = 0, pox = Poy = O.The point line is distribution for 
t = 0, the dashed line is distribution for t = 50T , no modulation e = 0, and the solid line is the 
distribution for t = 50T , e = 0.7. 
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FIG. 9. The probability distributions \ip{x,0)\'^ for y = 0. The initial wave function is the least 
uncertainty state, ax = cry = 0.1, xq = yo = 0, pox = Poy = 0.1. The point line is distribution for 
t = 0, the dashed line is distribution for t = 50T , no modulation e = 0, and the solid line is the 
distribution for t = 50r , e = 0.7. 
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